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Last year the E866-group of the Fermilab measured the xp dependence of J/\l/ and 
^ ! suppression in pA collisions. We discuss two of the effects found in that experiment with 
i/^ [ regard to color coherence effects: the different suppression of the J/\l/ and the "if' at xp < 
O I and the significant suppression of both at large xp- The small xp regions is dominated by 
Q ' fully formed charmonium states and thus enables us to discuss the formation time and the 
0^ ■ cross section of the different charmonium states. In the large xp region the interaction 
^ ' of the charmonium states with nuclear matter has to be described by partonic degrees of 
freedom, because in that kinematic domain the formation time is much larger than the 
nuclear radii. The understanding of this region will be crucial for the interpretation of 
^ ' the data of the future heavy ion colliders RHIC and LHC. 

c ■ 



1. Introduction 



From considerations on formation times the interaction of charmonium states with a 
\ nucleus in a pA collision can be roughly divided in three different kinematical regions. 
The formation time is given by |]T| 

1 2t> 

Here, Ey^/ j/^ are the energies (masses) of the and the J/\l/ and 7 is the 

Lorentz-factor. The formation time of a should be larger than that of a J/\l/ by the 
ratio of the radii of these states, which is ~ 2. In a recent letter the attempt to extract 
the formation time from the experimental data on e^e~ QQ has been made. The 
authors found tf{J/'^) = 0.44 fm/c-7, = 0.91 fm/c-7 which is slightly larger than 

the values used here. Using eq. (|]) in the rest frame of the nucleus leads to the following 
three kinematical regions: 

• hadronic: tf < 1.8 fm/c 

The formation time is smaller than the average inter-nucleon distance in the nucleus. 
Only fully formed charmonium states interact with the nucleons. In pA collisions 
at Eiab = 800 GeV this means xp < —0.3. 
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Table 1 

The non-perturbative charmonium nucleoli cross 
section (1) with the transverse size b calculated 
with the wave function from ||^ and (2) from [0]. 

• mixed: 1.8 fm/c < tf < 2 ■ Ra 

The formation time is larger than the average distance of the nucleons in the nuclei 
but smaller than the diameter of the nucleus. The charmonium states already can 
interact within their formation time but some of them are still formed in the nucleus. 
In pA {Eiab = 800 GeV) this means —0.3 < Xf < (the upper limit depends of 
course on A, x/ < is for Be). 

• partonic: 2 ■ Ra < tf 

All charmonium states are formed beyond the nucleus. Hence the interaction has 
to be treated on the parton level. 

In this work we discuss the hadronic and the mixed region. In section ^ we describe the 
model and present the results. A brief summary is given in section |^. 

2. Model and Results 

In we showed that the cross section of a fully formed J/\E' as predicted by pQCD |Q is 
much smaller than the cross section extracted from photoproduction processes at Eiab = 
20 GeV (here we have 7 = 6 and thus are in the hadronic region). We concluded 
that the cross section of charmonium states at Fermilab energies is dominated by a non- 
perturbative contribution. At the same time the pQCD contribution to the cross section 
increases rapidly with energy. It will be relevant at RHIC and may even dominate at 
LHC. 

To calculate the non-perturbative contribution to the charmonium nucleon cross section 
(shown in Tab. |I]) we used the parametrization a{b) = C ■ IP'. Here, h is the transverse 
diameter of the charmonium state (transverse to the beam direction) and C is a constant. 
This model is capable to describe the total cross sections of vr — A^, K — N, (p — N 
and pp collisions. We adjusted the constant C = a{n)/b'^{n), i.e. a{7i) = 25 mb and 
6^(7r) = 0.96 fm^ to the cross section and the size of the pion. (Note that the size of 
the pion is estimated from form factor measurements and thus is based on the idea of a 
non-relativistic continous charge density. Thus it may appear dangerous to use it as a 
spatial distribution of current quarks.) 

The transverse sizes of the charmonium states has been calculated with non-relativistic 
wave functions taken from The results of this calculations differ by ^ 25%. This 

difference illustrate the dependence of the prediction on the charmonium model. The 
subscripts of xim are the quantum numbers of the spherical harmonics, the angular mo- 
mentum dependent part of the wave functions. The x mesons have two different cross 
sections because the transverse size of P-states (/ = 1) depends on the third component 
of the angular momentum m: O"(xii)/o"(xio) = 2. This angular momentum dependent 
absorption leads to a polarization of the x-meson in pA and AB collisions. 

We assume in this work that the production of cc-pairs is a hard process, i.e. a{pA — »• 
cc) = A-a{pp —>■ cc). In the kinematic region discussed here, this assumption is consistent 
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with the data (for a recent review see e.g. P]) but may be inappropriate at large Xi;' 0. 
Additionally we assume that the suppression of charmonium states is only due to their 
final state interactions. In particular, we neglect the nuclear modifications of the parton 
distribution functions. 

In the following we calculate the survival probability of a state X in pA-collisions 
defined as 5*^ = a{pA X)/{A ■ a{pN — > X)). Thus S* = 1 for the production of cc- 
pairs (the total charm cross section), with the assumptions made above. Our assumptions 
allow to calculate the survival probability within the semiclassical approximation of the 
Gribov-Glauber model as in ||1CI|| : 



A-jJd^Bdz p{B, z) ■ exp [- axp{B, z')dz'^ , (2) 

where B is the impact parameter, z the beam direction, p the density of the nucleus and 
ax the cross section for a XN collision. 

Now, we can calculate 5* in the kinematical region that we called hadronic in the 
introduction. However, we also want to discuss the mixed region and therefore we have to 
describe the XN cross section within the formation time of the state X. We assume that 
the state X is produced as a small color singlet cc pair that expands during its formation 
time. Interesting physics related to the propagation of color octet QQ pairs is beyond our 
scope since it requires methods beyond the semiclassical approximation. For the region 
of negative xp the propagation of color octets should be suppressed because its energy is 
insufficient to separate color from the quark-gluon environment. The cross section then 
reads [O]: 



— z 

a{z,z')x = o"(-2) H {ax — <y{z)) for z' — z < tf 

a{z,z')x = CTx otherwise. (3) 

Here, a{z, z')x is the cross section of the state X produced at the position z after it 
moved to z'\ cr(z) = 1 mb is the cross section at the production point z and ax is the 
cross section of the state X after the formation time elapsed. That means a{z^ z')x oc t 
and thus b oc ^/t as it should be for the expansion of wave packages in nonrelativistic 
quantum mechanics This behaviour is called quantum diffusion. 

In Fig. H R{W/Be) := S{pW)/S{pBe) is plotted versus xp for the J/-^ and the 
Note that a large fraction of the J/\I''s comes from the decay of \E'' and therefore: 

^ = 0.6 ■ (0.92 ■ + 0.08 ■ 5*') + 0.4 ■ , (4) 

where S*"^/* is the genuine survival propability of the J/\£', while S is the survival prob- 
ability that takes into account the J/\E''s from the decay of the higher resonances. This 
calculation was done with the cross sections denoted (2) in Tab. |l| and the data from 



Ref. |12|. At < 0, we find reasonable agreement between data and calculation. At 



Xi? = we have Ic = 20.6 ■ 0.3 fm > 2 ■ Rbc and thus eq. (y) is not longer valid. 
3. Summary 

• The pQCD contribution into a{J /'^ — N) is too small at SPS-energies to explain 
a{'^A J/^)i and increases rapidly with energy - this could be important for 
RHIC and LHC. 
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Figure 1. 

The ratios R{W/Be) : = 
S{pW)/S{pBe) of the J/* and 
the are plotted versus xp at 
Eiah = 800 GeV. The used model 
(see text) is only valid for Xi;' < 0. 
In this region the calculation is in 
good agreement with the data. 
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• Motivated by the pQCD cross section, the nonperturbative model for (T{qq) = c-h"^ 
predicts four different types of charmonium nucleon cross sections. This reflects the 
fundamental property of QCD, where interactions depend strongly on the scale. 

• QCD predicts different cross sections for the Xim-niesons, because high energy in- 
teractions are dominated by the transverse size: — N) 2 • (y[xu) ~ N). 

• The production of small cc pairs and their expansion are relevant for the \[^':J/\1/- 
ratio, which is weakly dependend on A at larger Lorentz-factors. 

• The assumption that the production of cc-pairs is hard and not soft could be tested 
by measuring the ^-dependence of the total charm cross section. This could also 
provide insight in nuclear shadowing and parton energy loss at large Xp. 
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